Form Factors in D Meson Decays 



Dae Sung Hwang^ and Do- Won Kim'- 6 ) 
a: Department of Physics, Sejong University, Seoul 143-747, Korea 
b: Department of Physics, Kangnung National University, Kangnung 210-702, 

Korea 



Abstract 

We study the dT/dq 2 spectra and the branching fractions of the D meson exclusive 
semileptonic decays with the lepton mass effects into consideration. We investigate 
their sensitivity to form factor models, and find that the decays to a pseudoscalar 
meson and a lepton pair are sensitive to the property of the form factor Fi(q 2 ), and 
those to a vector meson and a lepton pair to the form factor Ai(q 2 ). We also an- 
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alyze the experimental results of the branching fractions B(D — > K~(or 
and B(D° — > K~(ot n~) e + z/), and show that it is implied that Fi(q 2 ) is of dipole 
type, instead of simple pole type which is commonly assumed. 
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1. Introduction 



The CP-violation phenomenon is discovered only in the Ki — > txtx decay and 
the charge asymmetry in the decay Ki — > 7r ± / =F z/ for more than 30 years. The 
B-factories at KEK and SLAC are under construction for the discovery of CP- 
violation in the B meson system. The mechanism of CP-violation through the 
complex phase of the CKM three family mixing matrix [|l| is presently considered 
standard for the CP-violation. In order to measure the CKM matrix elements 
accurately, it is important to know the hadronic form factors of the transition ma- 
trix elements reliably. For the heavy to heavy transitions the heavy quark effective 
theory provides good informations for the form factors. However, for the heavy 
to light transitions the understanding of the form factors is still limited and this 
fact hinders the extractions of the CKM matrix elements from experimental results 
significantly. At the same time, we will be able to have important clues for the 
internal structures of hadrons by knowing these form factors well. Semileptonic 
decay processes are good sources for the knowledge of the form factors both exper- 
imentally and theoretically, and the lepton mass effects in heavy meson exclusive 
semileptonic decays were studied by Korner and Schuler 0. 

We derive the formulas for dT/dq 2 with non-zero lepton mass in the forms 
which are efficient to study the form factor dependences. This formula for the 
pseudoscalar to pseudoscalar transition was also given by Khodjamirian et al. [|J. 
By using these formulas we study the dT/dq 2 spectra and branching fractions of 
the exclusive semileptonic D meson decays: D° — > K(*^e + v, D° — > 
D° — ► 7r _ (or p~) e + v and D° — > 7r _ (or p~) p + u. In this anslysis we employ three 
models of form factors and show how the results are influenced by the difference 
of form factors. For the decays to a pseudoscalar meson and a lepton pair, the 
results are sensitive to the property of the form factor F\(q 2 ). For the decays 
to a vector meson and a lepton pair, the results are sensitive to the form factor 
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Ai(q 2 ), and not to the other ones (A (q 2 ), A 2 (q 2 ) and V(q 2 )). We also analyze 
the experimental results of the branching fractions B(D° — ► K~(or 7r~) 7r + ) and 
B(D° — > i^~(or 7r~) e + i/), and show that it is implied that the form factor Fi(q 2 ) 
of the D to K and the P to 7r transitions are of dipole type, instead of simple pole 
type which is commonly assumed in the studies of the D meson decays. 

2. Semileptonic Decays of Heavy Mesons 

From Lorentz invariance one finds the decomposition of the hadronic matrix 
element for pseudoscalar to pseudoscalar meson transition in terms of hadronic 
form factors: 

< P(p)\J„\P(P) > 

= (P + p)J + (q 2 ) + (P-p)J_(q 2 ) 

f, s M 2 - m 2 \ „ , 9x M 2 -m 2 „ . 9x 
= ((P + p) M ^ ^)^i(? 2 ) + ~ 2 q»Fo(q 2 ), (1) 

where J M = ^^(l — 75) q- We use the following notations: M represents initial 
meson mass, m final meson mass, m; lepton mass, P initial meson momentum, p 
final meson momentum, and g M = (P — The form factors Fi(q 2 ) and F (q 2 ) 
correspond to 1~ and + exchanges, respectively. At q 2 = we have the constraint 
Pi(0) = Po(0), since the hadronic matrix element in (P is nonsingular at this 
kinematic point. 

The q 2 distribution of the semileptonic decay D° — > K~l + v is given in terms 
of the hadronic form factors Fi(q 2 ) and F (q 2 ) as: 

dT(D° — > K~l + v) Gl . l2 . 2w m,\ 2 

[ (K(f)f (1 + i^) |F,tf )| 2 + M 2 (1 - 1 ^F,!,')!'] , 
where K(q 2 ), momentum of the final meson in the D meson rest frame, is given by 
K(q 2 ) = ^({M 2 + m 2 -q 2 ) 2 - AM 2 m 2 ) * , (3) 
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and the physically allowed range of q 2 is given by 

m] <q 2 <(M - m) 2 . (4) 
For mi — 0, @ is reduced to the commonly used well-known formula: 

dT{D ° 7/~^ } = ^ IK/ (^(g 2 )) 3 |F l( g 2 )| 2 , (5) 

and < q 2 < (M — m) 2 . We note in (|5|) that only F\(q 2 ) contributes for mi = 0, 
however, for mi ^ F (q 2 ) also contributes as we can see in (0). 

From Lorentz invariance one finds the decomposition of the hadronic matrix 
element for pseudoscalar to vector meson transition in terms of hadronic form 
factors: 



<V(p)\J ll \P(P)> 



e-( P )((M + m) M >) - 2^^) + ^Mq 
+^e^^-V(q 2 )), (6) 



M + m 
where £0123 = 1 and 

2mA (q 2 ) = (M + m)A 1 (q 2 ) + g A 2 (q 2 ) + j^—A^q 2 ). (7) 

M + m M + m 

The form factors V(q 2 ), Ai(q 2 ), A2(q 2 ) and Ao(q 2 ) correspond to 1~, 1 + , 1 + and 
0~ exchanges, respectively. At q 2 = we have the constraint 2mA (0) = (M + 
m)Ai(0) — (M — m)A 2 (0), since the hadronic matrix element in (|16|) is nonsingular 
at this kinematic point. 

After a rather lengthy calculation, the q 2 distribution of the semileptonic decay 
D° — > K*~l + v is given in terms of the hadronic form factors Ai(q 2 ), A 2 (q 2 ), A 3 (q 2 ) 
and V(q 2 ) as [|: 

d T(D° ^ K-l+v) _ Gl l K(rt(1 _<, x (8) 



- 32tt 3 1 1 M 2 ^ J " q 2 



3 g z q z 2 



+Re(A 1 (g 2 K(g 2 ))[- M2 % * &MK)\\ - % + 2(MK) 2 ^ 

m z S q z q z 

1 



+ ^(M 2 + m 2 - q 2 )m 2 } + (M 2 - m 2 + q 2 )m 2 } 
+ l^(^)| 2 (M + 1 m)2m2 (^) 2 ^(^) 2 ( 1 <MK) 2 ^ + 2M 

+\V(Q 2 )\ 2 777^ , 2 [l(MK) 2 (l - 4) + 4(MiT) 2 4] 



2 m 2 l 



(M + m) 2L 3 v 7 v q 2 1 v ' g 2 

q 2 1 
( 

(M + m) 2 m 2 2 



+1^3(g 2 )l 2 777 ~ {MKfm 2 



-Re(A 3 (q 2 )A* 2 (q 2 )) 1 2 (M 2 - m 2 + q 2 )(MK) 2 m 2 

[M + m) z m z 



+Re(A 3 (q 2 )Al(q 2 ))^(MK) 2 m 2 }. 
m z 



When we take mi — * in (|8|), it agrees with the formula for = given in Refs. 

5 W (9) 

where 

i zLA/f 2 ^ 2 

H\q 2 ) = -4^((M 2 -m 2 -g 2 )(M + m)A(g 2 )-^^A 2 (g 2 )),(10) 

tf±(g 2 ) = -{(M + m)A 1 (q 2 )T^^V(q 2 )). 

In the case of the B to D meson (heavy to heavy) transition, the heavy quark 
effective theory (HQET) gives the useful relations between the relevant form factors 

B 

Fl (q 2 ) = V(q 2 ) = A (q 2 ) = A 2 (q 2 ) = F(y), (11) 

T-i / 2\ A ( 2x 2y/Mm y + 1 

m ) = am ) = — Hvi 

where y = (M 2 + m 2 — q 2 ) / [2Mm) = E D {*)/m (E D {*) is the energy of meson 
in the B meson rest frame), and F(y) is a form factor which becomes the Isgur- 
Wise function in the infinite heavy quark mass limit. When we use the relations 



flllD, for mi = the formula ([8]) becomes the well-known formula for the B to D* 
transition: 

dT{B° -> D*n~v) G 2 F 



dq 2 48vr 3 



l^ 2 m 3 (M - m) 2 ^/y 2 -l (y + I) 2 x 



where r = m/M. 

3. L»° K^-l+is 

For the form factors concerned with the exclusive semileptonic decays of D 
meson, we can not use the relations ([TT|) of the HQET. Therefore, in the study of 
D meson decays we use models for form factors. The pole-dominance idea suggests 
the following q 2 dependence of the form factors [§] : 

\ mi * 

where nj i and mj i are corresponding power and pole mass of the form factors fi(q 2 ), 
respectively. The WSB model |5J adopts Uf i = 1. However, the exact values of nf i 
are not known. The relations ( |TT| ) of the HQET gives the following approximate 
relation among the powers of the form factors for the heavy to heavy transitions: 

n Fl =n v = n Ao = n A2 = n Fo + l = n Al + l. (14) 

Non-perturbative analysis of QCD || suggests the same relation as (|T4|) for the 
form factors of the heavy to light transitions. The lattice calculations also show 
that the form factors F\, V and Aq are more rapidly increasing functions of q 2 
than the form factors F and A\ P, [T^] , which favors the relation (|14[) . Therefore, 
we will adopt two other models incorporating the relation (|14]) , as well as the 
WBS model which assumes nj. = 1 §, for the study of the exclusive semileptonic 
decays of D meson. We organize in Table 1 the values of the powers n/ s of the 



three models which we use in this work. For the values of the pole masses and 
those of the form factors at q 2 = 0, we use the values organized in Table 2 and 
3, which were given by Wirbel, Stech and Bauer ||. Their precise values are not 
known and they should be different for each of the three models. However, their 
exact values are not crucial for our work of clarifying the model dependences of 
the dT/dq 2 spectra and the branching fractions. 

For D° — > K~l + v, we use the formula (g) with non-zero lepton mass, instead of 
the commonly used formula (0) which is true for zero lepton mass. The obtained 
dT(D° — > K~l + v)/dq 2 spectrum and branching fractions are presented in Figure 

I and Table 4, for each of the three models we adopt in this work: WSB, Model I 
and Model II explained in Table 1. We find that the shape of spectrum of Model 

II is different from those of WSB and Model I in Figure 1. That is, the value 
of rip 1 determines the shape of spectrum. The experimental result of the E687 
Collaboration for this spectrum was given in Fig. 3 (a) of Ref . |13| , and the shape of 
their spectrum favors Model II better than WSB and Model I. In the experimental 
extraction of the value of F® K (0), the simple pole of the form factors has been 
commonly assumed [11, [l!|. Under this assumption, F® K (0) = 0.75 ± 0.02 ± 0.02 



was extracted JTT[| from the experimentally measured branching fraction B(D° — s 
K'e+iy) = (3.68±0.21)xl0~ 2 . (In Ref. pf, F° K (0) = 0.76±0.03 was presented.) 



However, if we assume in this extraction Model II which has the dipole form factor 
for Fi(q 2 ), we would get 



F? K (0) = 0.75 x — —- = 0.75 x 0.85 = 0.64 (15) 
1 v; V 4.78 x 10- 2 v ; 

for the mean value. In (|lf|) we used our results in Table 4 of the branching fraction 
B(D° — > K~e + v) for WSB (rip 1 = 1) and Model II (np 1 = 2), which were obtained 
by using the same value of F® K (0). That is, the experimentally extracted value 
of F® K (0) is much dependent on the form factor models used in the analysis. In 
reality, it is not yet established which type of form factors is the right one. 

For D° — > K*~l + v, we use the formula (ffi) with non-zero lepton mass. The 



obtained spectra and branching fractions are presented in Figure 2 and Table 5. 
We find that the results of WSB and Model II are almost the same, and they are 
significantly different from the results of Model I. This fact implies that the value 
of tia 1 mainly determines the spectra and branching fractions of D° — > K*~l + v. 

4. D° -> 7f-(or p~) l+u 

For D° — > n~l + v, we use the formula @ with the replacement of V cs by V C( i- 
The obtained spectra and branching fractions are presented in Fig. 3 and Table 

6. We find that the branching fractions of Model II in Fig. 3 are about twice 
those of WSB and Model I. Therefore, in case that we determine the value of 
F® w (p) from an experimentally measured branching fraction of D° — > tt~1 + v, the 
value of Ff n (0) determined with Model II will be about l/v2 times its value 
determined with WSB or Model I. From Table 4 and 6, we also find that the ratio 
B(D° — > 7i~l + v) /B(D° — ► K~l + v) from Model II is pretty bigger than that from 
WSB or Model I, but its present experimental result 0.101 ± 0.020 ± 0.003 |TJ can 
not discriminate them yet. 

For D° — > p~l + v, we use the formula (§) with the replacement of V cs by V c d- 
The obtained spectra and branching fractions are presented in Fig. 4 and Table 

7. We find in Fig. 4 and Table 7 that the results of WSB and Model II are almost 
the same, and they are much different from the results of Model I. Therefore, like 
the D° — > K*~l + v case, the property of the form factor A\(q 2 ) mainly determines 
the spectra and branching fractions, and the results are not sensitive to the other 
form factors (A (q 2 ), A 2 {q 2 ) and ^(g 2 )). 

5. Implications of Experimental Results 

In this section we compare the exclusive semileptonic decays and the two-body 
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hadronic decays. We start by recalling the relevant effective weak Hamiltonian for 
the two-body hadronic decay D° — ► K~tt + : 

H eS = ^V: s V ud [CMOi + C 2 (fi)0 2 ] + H.C., (16) 

1 = (uTPd)(sT p c), 2 = (sTPd)(uT p c), (17) 

where Gp is the Fermi coupling constant, V cs and V ud are corresponding Cabibbo- 
Kobayashi-Maskawa (CKM) matrix elements and V p = 7 p (l — 75). The Wilson 
coefficients C\(fi) and C 2 (n) incorporate the short- distance effects arising from 
the renormalization of 7i e R from /i = to /i = 0(m c ). By using the Fierz 
transformation under which V — A currents remain V — A currents, we get the 
following equivalent forms: 

1 

N c 
1 

N c 

where N c = 3 is the number of colors and T a 's are SU(3) color generators. The 
second terms in ([18]) involve color-octet currents. In the factorization assumption, 
these terms are neglected and 7i e s is rewritten in terms of "factorized hadron 
operators" 0: 

H cS = ^V*V ud (a 1 [uT'>d] H [sr p c] H + a 2 [sT p d} H [uT p c} H ) + H.C., (19) 

where the subscript H stands for hadronic implying that the Dirac bilinears inside 
the brackets be treated as interpolating fields for the mesons and no further Fierz- 
reordering need be done. The phenomenological parameters a\ and a 2 are related 
to Ci and C 2 by ai — C\ + jfC-z and a 2 = C 2 + ^-Ci- The numerical values of ai 
and a 2 for D meson decays are given by Jl5| 



C 1 1 + C 2 2 = (Ct + — C 2 )0 1 + C 2 {sTPT a d){uT p T a c) 
(C 2 + -^COOs + C^uYP^d^Y^c) 



a x = 1.10 ±0.05, a 2 = -0.49 ±0.04. (20) 

For the two body decay, in the rest frame of initial meson the differential decay 
rate is given by 
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[(M 2 -(m 1 +m 2 ) 2 )(M 2 -(m 1 -m 2 ) 2 )]^ 
IPil 2M ' (22) 

where M is the initial meson mass, mi and m 2 the final meson masses, and pi the 

momentum of one of the final mesons in the initial meson rest frame. By using 

flip, (|T9| ) and < 1 1 7r~ (g) >= iq^f n -, ( pi]) gives the following formula for the 

branching ratio of the process D° — > K~n + : 



R, n" ?■,---+- -- l ^F m2 D \2 | T/ |2 1 m D 2 /J | T/ rpDK ( 2\ |2 



On the other hand, from (|5p and (|13|) the branching ratio — > K~e + z/) is 



given by 



B(D°^K-e+u) = (^?^^\V cs F°«(0)\*xI™, (24) 



where the dimensionless integral I DK is given by 



((I + ^ -r)*- I-" 

= / 5- 

/n m„ 

-X 



/ ^ 2 N2raf ,'" g/ ■ (25) 



f 1 _ m D 

In the above, we neglected the electron mass. From ( |2"3"P and (^) we have 
8(£>° - AT-e+i/) ~ 1 m 2 , V 1 





to 






m 2 D ' 








-) 2 ) 





a? 



= 0.225 x— L:, (26) 
where we used the fact F® K (m 2 ) ~ ^^^(0) and the following experimetal values 



[0: m D = m D o = 1.8646 ± 0.0005 GeV, m K = m K - = 493.677 ± 0.013 MeV, 
m w = m w+ = 139.56995 ± 0.00035 MeV, f w = f w+ = 131.74 ± 0.15 MeV and 
V ud = 0.9753 ± 0.0008. 
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When we use the experimental results B(D° — > fT 7T + ) = 3.85 ± 0.09 % and 
£(£>° -> i^-e+z/) = 3.66 ± 0.18 %, © gives 

/^[Expt.] = 0.213 (1 ±0.054) aj = 0.258 (1 ± 0.054) (1 ± 0.091) 

= 0.258 (1 ±0.106) = 0.231 - 0.286, (27) 

where we used the value of a± given in (|20"|). On the other hand, when we calculate 
jdk from (|25|) with = 2.11 GeV geven in Table 2, we obtain the 

following results: 

I DK [n Fl = 1] = 0.195 for n Fl = 1, 

7 Dir [n Fl = 2] = 0.267 for n Fl = 2. (28) 

From (^) and fl28|), we find that the experimental results of B(D° — > K~tt + ) and 
-> K'e+v) imply n Fl = 2. 
In the same way as the above, for the D to 7i transition we get the formula 

B(£>° - Tr'Tr*) r 2|T/ |2 £ , m^2 



— — — _ = 6tt \V ud \ — 1 1 

£>(.D U — > 7T e+z/) m|p v 



V ud F n DK (0)\ 2 I D * 



m D ' | V ud Pq 
= 0.275 x , (29) 

where we used the fact F® 7T (m 2 r ) ~ F® n (0), and the dimensionless integral I Dn is 
given by 

r (i-**-r (a + # - x) 2 - 4^) 

2n Fl 



3 

2 _ X )2 _ 4 m2 r ^ 5 

/ T ' ^ ■ (30) 



o 



When we use the experimental results B(D° — > n 7T + ) = (1.53 ± 0.09) x 10 3 
and B(D° -> 7T- e + z/) = (3.7 ± 0.6) x 10~ 3 , (p§ gives 

J D ^[Expt.] = 0.665 (1 ± 0.173) a\ = 0.804 (1 ± 0.173) (1 ± 0.091) 

= 0.804 (1 ±0.195) = 0.648 - 0.961, (31) 
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where we used again the value of a% given in (|20|). When we calculate I D7T directly 
from ( p^D with m Fl = 2.01 GeV geven in Table 2, we obtain the following results: 

I D7T [n Fl = 1] = 0.385 for n Fl = 1, 

I D7T [n Fl =2] = 0.783 for n Fl = 2. (32) 

The results in (|32| ) show that the value of 7 D?r is more sensitive to n Fl than that of 
I DK . From fl3~l|) and fl32|), we find that the experimental results of B(D° — > 7r~7r + ) 
and £>(-D° — > 7r _ e + z/) imply = 2. 



6. Conclusion 



We studied the D meson exclusive semileptonic decays with the lepton mass 
effects into consideration, and investigated their sensitivity to form factor models. 
The results show that the decays to a pseudoscalar meson and a lepton pair are 
sensitive to the property of the form factor Fi(q 2 ), and those to a vector meson 
and a lepton pair are sensitive to the form factor Ai(q 2 ) and not to the other ones 
(A (q 2 ), A2(q 2 ) and V(q 2 )). Concerned with the experimental extraction of the 
form factor values at q 2 = 0, Ff K (0) = 0.75 ± 0.02 ± 0.02 has been obtained 
from the experimental result B(D° — > K~e + v) = (3.68 ±0.21) x 10~ 2 by assuming 
n Fl = 1 [|11[|, however, its mean value is modified to 0.64 if n Fl = 2 is adopted. 
In reality, the value of n Fl is not known presently. The experimental extraction of 
the reliable value of Ff K (0) is very important for the determinations of the fun- 
damental parameters, and also for the reason that the experimentally extracted 
value of Ff K (0) is compared with its values obtained by different theoretical cal- 
culations such as quark models, lattice QCD, and QCD sum rules [|ll], [l!|. We 



emphasized that this extraction is very much dependent on the value of n Fl . In 
this context, we also analyzed the experimental results of the branching fractions 
B(D° -> K-(oi 7T-) 7T+) and B(D° if- (or 7r~) e + u), and showed that it is 
implied that the form factor Fi(q 2 ) of the D to K and that of the D to ir tran- 
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sitions are of dipole type {np 1 = 2), instead of simple pole type (uf 1 = 1) which 
is commonly assumed in the studies of the D meson decays. We think that care- 
ful studies on the form factors are important not only for the extractions of the 
CKM matrix elements but also for the understanding of the internal structures of 
hadrons. 
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n Fl 


n Fo 


n v 


n Al 


n A 2 


KA 


WSB(pole/pole) 


1 


1 


1 


1 


1 


1 


Model I(pole/const.) 


1 





1 





1 


1 


Model II(dipole/pole) 


2 


1 


2 


1 


2 


2 



Table 1: The values of the power of pole for the three models used in this paper. 



Current 


to((T) 


m(l-) 


m(0+) 


m(l+) 


Relevant Form Factors 


A 


Fx V 


i^o 




sc 


1.97 


2.11 


2.60 


2.53 


dc 


1.87 


2.01 


2.47 


2.42 



Table 2: The values of pole masses (GeV) used in numerical calculations. 





Fi(0) = F (0) 


V(0) 


Mo) 


A 2 (0) 


A (0) 


D -> ir (,,t) 
L> ^ tt(p) 


0.762 
0.692 


1.226 
1.225 


0.880 
0.775 


1.147 
0.923 


0.733 
0.669 



Table 3: The values of the form factors at q 2 = used in numerical calculations. 
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B(D° -> K-e+v) 


B(D° -> 


B(e+) : B(^) 


WSB 


3.49 x 10" 2 


3.41 x 10" 2 


1 : 0.98 


Model I 


3.49 x lO- 2 


3.38 x 1(T 2 


1 : 0.97 


Model II 


4.78 x 10" 2 


4.67 x 10" 2 


1 : 0.98 



Table 4: The obtained branching fractions and their ratios for D° — > K 





B(D° -> if— 


e+v) 


B(D° -> 




s( e +; 


: 


WSB 


3.88 x 10" 


-2 


3.68 x 10 


-2 


1 


0.95 


Model I 


3.28 x 10" 


-2 


3.10 x 10 


-2 


1 


0.95 


Model II 


3.85 x 10" 


-2 


3.66 x 10 


-2 


1 


0.95 



Table 5: The obtained branching fractions and their ratios for D° — > K* 





B(D° -> vre+z/) 




B(e + ) : B(ji+) 


WSB 


2.92 x 10" 3 


2.88 x 10" 3 


1 : 0.99 


Model I 


2.92 x lO- 3 


2.86 x 10~ 3 


1 : 0.98 


Model II 


5.94 x 10" 3 


5.86 x 10" 3 


1 : 0.99 



Table 6: The obtained branching fractions and their ratios for D° — > n l + v. 





B(D° -> pe+z/) 


B(D° -> pp+z/) 


B(e+) : £(p+) 


WSB 


2.77 x 10" 3 


2.66 x 10" 3 


1 : 0.96 


Model I 


2.19 x 10~ 3 


2.10 x 10~ 3 


1 : 0.96 


Model II 


2.77 x 10" 3 


2.66 x 10" 3 


1 : 0.96 



Table 7: The obtained branching fractions and their ratios for D° — > p Z + z/. 
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Figure Captions 



Fig. 1. (l/T tot )(dT/dq 2 ) of (a) D° -> K'e+u and (b) D° -> fT>+z/ for three 
models: solid line for WSB, dashed for Model I, and dotted for Model II. Solid 
and dashed lines almost overlap. 

Fig. 2. {l/T tot )(dT/dq 2 ) of D° -> K*- e + u and £>° -> K*-fi + u for three models: 
solid line for WSB, dashed for Model I, and dotted for Model II. Solid and dotted 
lines almost overlap. 

Fig. 3. (l/T tot )(dT/dq 2 ) of (a) D° -> vr-e+z/ and (b) D° -> tt-/x+z/ for three 
models: solid line for WSB, dashed for Model I, and dotted for Model II. Solid 
and dashed lines almost overlap. 

Fig. 4. (l/Y tot )(dT/dq 2 ) of D° -> p~e+z/ and D° -> p>+z/ for three models: solid 
line for WSB, dashed for Model I, and dotted for Model II. Solid and dotted lines 
almost overlap. 
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